A Multidimensional Analogue of the 
Simpson 's Formula of Integral 



Kazuyuki FUJII * 
Department of Mathematical Sciences 
Yokohama City University 
Yokohama, 236-0027 
Japan 

Abstract 

The Simpson's formula is obtained by approximating the integral of a function on 
some interval by the integral of the quadratic polynomial determined by the function. 
However, a multidimensional analogue of the formula has not been given as far as we 
know. 

In this paper such a formula is given. Our formula is simple and beautiful, so it 
may be convenient in Mathematics or Mathematical Physics. 

1 Introduction 

Let y = f(x) be a continuous function defined on the interval [a, £>]. We want to calculate 
the integral f(x)dx. However, it is not easy to calculate the integral explicitly, so to make 



"E-mail address : fujii@yokohama-cu.ac.jp 



1 



; 



some approximation is a realistic way. We recall the Simpson's formula \j , see for example 
We divide [a, b] into equal 2n subintervals and set A = , and a = xq, b = %2n- 




a 
II 



X2n 



Now we study the integral on small interval [x2k-2, X2k] in the following. 




X2k-2 %2k-l X 2 k 

II II 

x 2 k-i - A X 2 k-l + A 



For 



we set 



V2k-2 = f{x2k-2), V2k-1 = f(%2k-l), U2k = f(^2k) 

A = (x 2 k-2, V2k-2) = (X 2 k-1 ~ A, V2k-2), 
C = (X 2 k-1, V2k-l), 

B = (x 2k , V2k) = {x 2k -i + A, y 2k ) 



1 The author uses the terminology Simpson's formula in place of Simpson's rule because it is nothing 
but an approximate formula (or method) 



for simplicity. 

A quadratic polynomial 



y = px + qx + r 



passing through three points A, B, C is determined uniquely. Namely, we have only to solve 
the simultaneous equations 

(x 2k -i - A) 2 p+ (x 2k -i ~ A) q + r = y 2k _ 2 , 
X \k-\V + x 2 k-\q + r = y 2 k-i, 

\2 



and the result is 



V 



{x2k-\ + A) p + (x 2k -i + A)q + r = y 2k 



y2k - 2|/ 2 fc-l + V2k-2 

2A 2 



V2k~V2k-2 ?/2fc - 2y2fc-l + V2k-2 



2A 

V2k-1 - 



y2k - V2k-2 

2A 



2A 2 

%2k-l + 



y2k - 2y 2 fe-i + y2 k -2 2 



2A 2 



x 



2k-V 



Therefore the function y — f(x) on the interval [x 2k - 2 , x 2k ] is approximated by y = px 2 + 
qx + r like 

B 



2/2fc-2 



1 


/ : 
















V2k-\ 



V2k 



%2k-2 %2k-l X 2k 

II II 

x 2k ^ x - A x 2k -i + A 



Then an important formula is obtained 



/•x 2 fc_i+A ^ 

/ {px 2 + qx + r)dx = — {y 2k - 2 + ^y 2k -i + Jfefc) 
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because 

LHS = ^{{x 2k _ l + Af-{x 2k _ l -Af} + l{{x 2k _ l + A) 2 -{x 2k _ l -A) 2 }+2rA 
= - P A 3 + 2( P x 2 2k _ 1 + qx 2k _ 1 + r)A 

= -pA 3 + 2y 2k _ 1 A 

A / 2 y 2k ~ 2y 2 fe-i + V2k-2 . a 
= y|2A x — + 6y 2k _ 1 

= ^(V2k-2 + %2fc-l + V2k) = RHS. 

From this we have the famous Simpson's formula 



/ f(x)dx = ^2 f{x)dx = ^2 f{x)dx 

Jo, k=1 J x 2 k-2 k= i Jx 2k _ 1 -A 

/ (px 2 + qx + r) dx 



k=i ux ^k-i- 

n 



= ^2^(y2k-2 + 4y 2 k-i + y2k) 

k=l 

b — a 



6n 



{yo + 4(2/i + y 3 H h Z/2n-i) + 2(y 2 + Vi H h 2/2^-2) + 2/2«} 



because A = V^- 

2n 



2 Generalized Simpson's Formula 

Let D be an n-dimensional cuboid 

D = [ai, h] x [a 2 , b 2 ] x • • • [a n , 6„] C R" 

and / be a continuous function on £>. What we do is to calculate the multidimensional 
integral 

/ I I f( Xl ' X2 '"' ^ x n)dx\dx 2 - ■ -dx n . (1) 
However, it is almost difficult, so we must satisfy only by approximating it as shown in the 
introduction. For that purpose we rewrite D as follows. 

D = [ax — Ai, a x + A x ] x [a 2 - A 2 , a 2 + A 2 ] x • • • x [a n - A n , a n + A n ] (2) 



where otj = (bj + aj)/2 and Aj = (bj — %)/2. See the following figure (n = 2). 



(a ± ,a 2 + A 2 ) 

(cti - Ax, a 2 + A 2 ) ; («i + Ai, a 2 + A 2 ) 



(ai - Ai, a 2 ) 



(«i - Ai,o!2 - A 2 




(a 1? a 2 - A 2 ) 



(ai + Ai,a 2 ) 



+ Ai,a 2 - A 2 ) 



Next we approximate / by a quadratic polynomial 



F(x ± ,x 2 , ■■■ ,x n ) = 

which is a natural extension in case of n = 1. 
For ji, j 2 , • • • ,jn e {-1,0, 1} we set 



it; 



= /(ai + JiAi, a 2 + j 2 A 2 , •••,«„+ j n A n ) 
= F(a x + jiAx, a 2 + j 2 A 2 , •••,«„+ j n A n ). 

It is of course 

j = 3" = |j{aui 2 ...i n }, 

see the figure above. Then we have 



(3) 



(4) 



Formula (Conjecture) 




/ F(xi,x 2 , • • • , x n )dxidx 2 ■ ■ ■ dx n 
Jd 



A!A 2 ---A ra 

3 n 



^ ^ ^ttO'l ,32,— Jn>0 



31,32,- JnG{-l,0,l} 

where j 2 , • • • , j n }o is the number of in {j 1 ,j 2 , ••' ,j n }- 



(5) 



For the case of n = 2 and n = 3 the formula is proved in the following. However, the 
proof of the general case is left to readers. 



From the formula we obtain a good approximation 

f (x\^ x 2 , ' ' ' ; x n }dx\dx 2 ' ' ' dx n 




D 



371 " 3x3: 

J1J2,- J«6{-1,0,1} 



if D is small enough. 



2.1 Proof for n = 2 



Let us prove (jSJ) for n = 2. From ([3]) we set for simplicity 



a 20 £ 2 + a 10 x + a o + (a 21 x 2 + sui + a i)y + (a 22 x 2 + a 12 x + a 02 )y 2 . 



See the figure once more. 



(at - Ai, a 2 + A 2 



(ai - Ai, a 2 ) 



(ai - Ai, a 2 - A 



(ai, a 2 + A 2 




(ai + A x , a 2 + A 2 ) 



(ai + Ai,a 2 ) 



(ai + Ai, a 2 - A 2 



(ai,a 2 - A 2 



6 



Then 



I / wdxdy 

ai-Ai Ja 2 -A 2 

Qi + Ai) 3 - (ai - Ai) 3 , _ 

^20 g h »io ^ 1- aoo-2Ai } 2A 2 

(a 1 + A 1 ) 3 -(g 1 -A 1 ) 3 ,. . L , , 
W + aii + a Qr 2Ai^x 



(a 2 + A 2 ) 2 - (« 2 - A 



2) 



+ S a22 



2 

(ai+Aa) 3 -(ai-Ai) 3 



(ai + A x ; 


) 2 - 


(ai 


-AO 2 


(ai + Ai; 


o 
Z 

) 2 - 


(ai 


-AO 2 




2 






(ai + A x ; 


) 2 - 


("i 


-AO 2 


2 



3 

(a 2 + A 2 ) 3 -(a 2 -A 2 ) 3 



= 4AiA 2 {a 20 o; 2 + ai «i + a 00 + (a 2 ia 2 + a u ai + a 01 )a 2 + (a. 22 a 2 + a^ai + a 02 )a 2 } 

+ 4A ^ 2 (fl2oA 2 + a 2i A 2 a 2 + a 22 A 2 a 2 + a 02 A2 + a^aiA^ + a 22 « 2 A 2 ) 

+ ^a 22 A 2 A 2 . (8) 

From (jlj) we have 9(= 3 2 )-data 

=w(ai + k Ai ,0:2 + ^2) (9) 

for k, I G { — 1, 0, 1}, so we want to express the integral above in terms of these data {wki}- 
Since 

«2o«i + «-io«i + «oo + («2i«i + aii^i + a Q i)a 2 + («22«i + ai2«i + «02)« 2 = ^o,o (10) 

is clear we must treat the remaining ones 

a 20 A 2 + a 2i A 2 a 2 + a 22 A 2 a 2 + a 02 A 2 + ai 2 «iA 2 + a 22 a 2 A 2 , 
a 22 A 2 A 2 . 

It is not difficult to show 

a 20 A 2 + a 2 iA(a 2 + a^A^ 



2 . \2., , X 2..2 W_ 1;0 - 2w 0i Q + W h0 

2 

A 2 , A 2 . 2 A 2 Wq,_i - 2Wq,q + Wp,i 

a 02 A 2 + ai 2 «iA 2 + a 22 a:iA 2 = . 
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so we have 



a 20 Ai + a 21 Aja 2 + a 22 A\a\ + a 02 Al + a 12 a x Al + a 22 a\A\ 
_ W-1,0 + Wi,o + wp,-i + wp.i ~ 4w ,o ^ 

Similarly, 



W-i-i + 10-1,1 + Wi _i + lOi. 



4 

= a20«i + fl io«i + «oo + («2i«i + anai + aoi)«2 + (022^1 + ai2«i + 002)0! 
+ a 2 oA^ + a 21 Afa 2 + a 22 A\al + a 02 Aj + a^aiA^ + a 22 a\A\ 

+ a 22 A^ 

w_i.o + Wl.O + u>o,-i + Wo,i - 4^0,0 . a2a 2 
= w ,o H ^ h a 22 A 1 A 2 



so we have 

Ct A^ A^ L '~ 1 ' — 1 l " I > _1 1 "'■M "V"'— i,u 1 "'ijU 1 "'u,- x 1 " / u,i/ 1 il ^u,u (12) 



2 2 _ iy_i_i + W-1,1 + wi _i + 101,1 - 2(w_i i0 + wi,o + + w o,i) + 4wo,o 



Therefore, substituting (TIU]). (fTTl) and (fl2]l into (ED we obtain 

-ai+Ai /•« 2 +A 2 



wdxdy 

Oil— Al J Ct2— A2 

2AiA 2/ 

= 4AiA 2 w ,o H g — (w-1,0 + wi.o + w -1 + io 0) i - 4u7 ,o) 

+ q 2 (w-i,-i + + + - 2(w_i j0 + Wi,o + ^o-i + Wo,i) + 4w ,o) 

AiA 2 r , -, 

= — - — + w-i,i + ^l-i + u>i,i + 4(10-1,0 + w 1)0 + w _i + w ,i) + 16w ,o} • 

y 

(13) 

This formula is interesting enough. 
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2.2 Proof for n = 3 

Let us prove (E]) for n = 3. From (jHJ) we set 

2 

w = y, z) = a ijk x l y 3 z k 

i,j,k=0 

= a 200 x 2 + a 10Q x + a 000 + (a 2 io^ 2 + «no^ + «oio)y + (0220a; 2 + ^120^ + a 020 )y 2 

+ {a 2 oiX 2 + a 101 x + a 00 i + (a 2 nx 2 + a m x + a ou )y + (a 221 x 2 + a 121 x + a 021 )y 2 } z 

+ {0202a; 2 + a 102 x + a 02 + (0212a; 2 + a U2 x + a 012 )y + [a 222 x 2 + 0122a; + a 022 )y 2 } z 2 . 

(14) 



Then 



CI1+A1 pa 2 +A 2 PO13+A3 

/ / wdxdydz 

a\— Ai J a 2 — A2 J as— A3 

pax+Ai ra 2 +A 2 2 1 

</ai-Ai Ja 2 -A 2 ij = o J 

Q1+Al , Q2+ a 2 ^ .| (g 3 + A3) 2 - (g a ~ A3) 2 



+ < / / 5^ a iii x ^ 

[Ja 1 -A 1 Ja 2 -A 2 j J= o 

j ,a 1+Al p + A 2 2 I + A 3 _ _ Ag 



oci— Ai Ja 2 -A 2 j J= o 



By (|SJ) some calculation gives 



9 



/■ai+Ai pa2+A2 raz+A-s 

/ / / wdxdydz 

j«i-Ai Jao—A? Ja-K—l 



' a\— Ai ^Q2~A2 Jag— A3 

= [4AiA 2 {a 2 ootti + aiooai + a oo + (0210^1 + a no«i + aoio)«2 + (0220^1 + a i2o«i + ao2o)a 2 } 

+ 4A ^ 2 (Q2ooAi + a 2 ioAia 2 + a^oA^ + a 2oA 2 + a 12 oot\A\ + 0220^1 Ag) 

4AiA 2 A 9 A 9l A 
+ — ^a 2 2oA?A2]2A 3 

+ [4AiA 2 {a 20 io;i + aioiai + a oi + {a 2 ua\ + o m ai + a n)o; 2 + {a 22 \a\ + flmai + a 02 i)a 2 } 

+ 4A ^ 2 (a 2 oi Ai + a 2 nAi« 2 + a 221 Ajal + amiA\ + o^iaiAj + a^ia^A^) 

4AiA 2 A 9 A 9l 
+ — jp^a 22 iA2A2]2a 3 A 3 

+ [4AiA 2 {a 202 a;i + a W2 ai + a 02 + (a 2 i 2 a: 2 + aii2«i + aoi2)«2 + (0222a 2 + ai 22 ai + ao22)a 2 } 

+ 1^1^ (a 2 o 2 A^ + a 2 i 2 A\a 2 + 0,222 A\a\ + a 022A 2 + ai 22 aiA2 + a 222 a\A\) 
o 

+ 4 ^a 222 A?A^2^A 3 + ^A^) 
= 8A!A 2 A 3 x 

[ {a 2 oo"i + aiooai + aooo + (a 2 io«? + anoai + a 010 )a 2 + (0220a 2 + a i2oai + 0020)^2} 
+{a 20 i«i + aioi«i + a oi + («2ii«i + «niai + a n)a; 2 + (a 221 a\ + a 121 a! + a 2i)a; 2 }a; 3 

+{«202a 2 + aio2«i + 0002 + (0212a 2 . + a ii2«i + 0012)0:2 + (0222a 2 + ai 22 «i + ao22)a 2 }« 2 ] 
8AiA 2 A 3 

[ {a 20 oAi + a 2 i Ai« 2 + a 220 A^2 + «02oA 2 + auo^Aj + a 220 a; 2 A 2 ! } 
+{a 20 iA^ + a 2 uA\a 2 + a 221 Ajaj + CI021A2 + a 121 axAl + a^ia^A^a^ 
+{a 202 A^ + a 2X2 A\a 2 + a 222 A\a\ + + a Y22 axA\ + a 222 a\A 2 2 }o? z \ 

8^]^^2^3 / a 2 a 2 a 2 a 2 a2a2 2\ 

+ g (a 220 A x A 2 + a 22i A 1 A 2 a 3 + a 222 A 1 A 2 a 3 ) 

8A!A 2 A 3 

{a 202 al + aio2«i + ^002 + (a 2 i2«? + flii2«i + a i 2 )« 2 + (a 222 aj + ai 22 «i + a 022 )a 2 }A3 
+ 8 ^ 1 '^ 2 ^ 3 ( Q 202 A\ + a 212 A\a 2 + a 222 AiO; 2 + a 022A 2 + ai 22 aiA2 + a 222 a:i A2)A 3 
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By rearranging terms we have 

r-ai+Ai /-02+A2 /-a3+A3 



wdxdydz 

Ot\— Ai Jot?.— A2 ^«3— A3 

= 8A!A 2 A 3 x 

[ {fl2oo«? + aioo«i + aooo + (a 2 io«i + a iio«i + aoio)a2 + (0220"? + ai20«i + «02o)«2} 
+{«-2oi«i + «ioi"i + 0001 + (a2ii«i + amai + a n)a 2 + (a 2 2i«i + ai2i«i + ao2i)« 2 }«3 

+ {a 202 «i + ai02«l + O-002 + (^212"! + GHl2ttl + «012)"2 + (a-222ttl + 0122^1 + O-022)a 2 }a 3 ] 

8A!A 2 A 3 

+ x 

3 

[ {a 20 o + ct2ioa2 + ^220^2 + ( a 2oi + 0-21102 + o 22 ia 2 )ct3 + (a 202 + o 2i2 a 2 + a 222 a 2 2 )al}A\ 
+{0020 + ai2oai + 0220"! + ( a 02i + flmQti + a22itti)a 3 + (a 02 2 + o i22 ai + o 222 ai)a 3 }A2 



3 J 



+{«oo2 + «io2«i + «202«i + (0-012 + aii2«i + a 2Y2 a^)a 2 + (a 02 2 + «i22«i + a 222 a 1 )a 2 }A 

8A!A 2 A 3 

+ x 

9 

{a 2 ^A\A\ + a 221 A\A\a z + a 222 A\A\a\ + a 202 A?A^ + a 2 i 2 A?a 2 A^ + a^A^A^ 
+a 02 2A 2 A3 + a X22 a\A\A\ + a^a^A^A^} 

a 222 AlA 2 2 Al (15) 



8A!A 2 A 3 a2a2a2 



27 

From (jU) we have 27(= 3 3 )-data 

^,m,n = w(ax + /Ai, a 2 + mA 2 , a 3 + nA 3 ) (16) 

for I, m, n G { — 1, 0, 1}, so we must express the integral above in terms of these data {w^ myn }. 
Since 

{a 200 ai + aioo«i + o 00 o + (a 210 a\ + ano^i + a 010 )a 2 + (a 220 a\ + a 120 a 1 + a Q20 )a 2 } 
+{0201^1 + Oioitti + a ooi + (o2ii«i + ctiiitti + 0011)^2 + (o-22iai + O121CH + a-02i)a 2 }o; 3 
+{a202«i + aio2«i + CJ002 + (a2i2«i + aii2«i + o i 2 )a 2 + {a 222 a\ + ai22«i + o 022 )a 2 }a 3 
= w ,o,o (17) 

we treat the remaining ones. Here we list formulas necessary to prove the main formula. 
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{a 20 + 0210^2 + fl220«2 + ( fl 201 + a211«2 + a 22 ia 2 )«3 + (a 2 02 + «212«2 + 0222a 2 )«3}Al 

1/ 

= ^(w-i.o.o + w i,o,o) - ^0,0,0, 

{a 20 + «120«1 + «220«1 + («021 + ai21«l + a 2 2iai)a3 + (a 02 2 + «122«1 + a222Q ! l)«3}A 2 
= ^(^o -1,0 + Wo,l,o) - ^0,0,0, 

{^002 + O-102«l + «202Ol + (^012 + Oll2«l + a^O^O^ + («022 + 0122^1 + «-2220l) a 2}^3 

= -(it;o,o,-i + wo,o,i) - wo,o,o (18) 



and 



and 



a 22 o^?^2 + 022iAiA|a 3 + a^A^A^ag + a^A^Ag + a^A^A^ + a^A^a^Ag 

+CI022A2A3 + a^aiA^Ag + a^a^A^A^ 
1, 

= -(w-i -1,0 + W-1,1,0 + wi ,-i,o + ^1,1,0 + W-1,0,-1 + W-1,0,1 + ^1,0,-1 + ^1,0,1 
+W _i _i + W + U>0,1,-1 + ^0,1,1) 

-(w-1,0,0 + u>i, ,o + ^0 ,-1,0 + ^0,1,0 + ^0,0,-1 + Wo,o,i) + 3wo,o,o (19) 



!/ 

= g(w-i -1,-1 + wi-1,-1 + + + + ^1 -1,1 + w i,i,i + 

1, 

--(W_l _l j0 + W-1,1,0 + ^1 -1,0 + Wl,l,0 + W-1,0,-1 + W-1,0,1 + Wl.0,-1 + ^1,0,1 
+Wq _i _i + W -1,1 + M>0,1,-1 + ^0,l,l) 

+ 2-( w ~i,o,o + ^1,0,0 + w -1,0 + w oxo + iwo,o-i + ^0,0,1) - ^0,0,0- (20) 

The proofs are long but straightforward, so they are left to readers. 
Therefore, substituting ifTTj) ~ (1201) into (IT51) we obtain 
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•ai+Ai /-Q2+A2 ras+A3 

wdxdydz 



J a\— Ai J«2- A2 Jots— A3 

= 8A1A2A3 u; ,o,o 



;AiA 2 A 3 f 1 1 
< -(lW-1,0,0 + ^1,0,0 + w -1,0 + IWO.1,0 + w , ,-i + w ,o,i) - 3w ,o,o r- 

8A!A 2 A3 f 1 , 

+ g S ^(lW-1,-1,0 + W-1,1,0 + ^1,-1,0 + ™1,1,0 + W-1,0,-1 + w-1,0,1 + w i,o,-i + w i,o,i 

+W + Wo -1,1 + W ,l,-1 + W ,l,i) 

-(w_i A o + Wi,o,o + w Q -1,0 + w ,i,q + Wofl-i + w , ,i) + 3u; ,o,o} 
SAxAsAg f 1 . 

+ 27 1 g ( u '-i -1 -1 + -1,-1 + ^1,1,-1 + w -i,i -1 + w -i -1,1 + wi -1,1 + + w_ 

1, 

--(10-1-1,0 + W-1,1,0 + Wi _i,o + Wl,l,Q + W-1,0 _i + W-i,o,l + Wi )0 ,-l + lWi,o,l 

+W ,-1,-1 + W -1,1 + W ,l,-l + w ,i,i) 
+ 2-(W-i,o,0 + Wi A o + Wo-1,0 + W Q)1 fl + W , ,-l + w , ,i) - w , , j 



AiA 2 A 3 



{w_i _i_i + Wx _i_i + W lt i _l + W_i,i _i + W-1,-1,1 + Wi -1,1 + + W_i,i,i 

+4(w_i _i )0 + W-1,1,0 + Wl -1,0 + Wi,i, + W-1,0,-1 + W-1,0,1 + Wi, -1 + Wi, ,l 

+W _i _i + W -1,1 + W ,i -1 + W ,i,i) 
+ 16O_i,0,0 + Wi, ,o + W ,_l,o + W ,1,0 + w , ,-i + W ,o,l) 

+64w ,o,o} • 



The proof is a bit complicated. 



2.3 General Case 



Let us consider the general case. Unfortunately, we have no method to calculate the general 
case at the present time. However, from the results (TlBl and (121j) we can conjecture the 
general formula as follows. 
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The formulas ({TBI and f|2Tjl can be rewritten as 

"«1+Al /-«2+A2 



wdxdy 

Ql- Ai J Ct2— A2 

= + w-i,i + wi,-i + wi,i + 4(w_i,o + ^1,0 + w 0) _i + w ,i) + 4 2 w ,o} 

and 

«q:i+Ai ca2+A2 /-a3+A3 

/ wdxdydz 

ai-Ai J a-2— A2 Ja?,— A3 

AiA 2 A 3 

+ + + + W-l -1,1 + + ^1,1,1 + ^-1,1,1 

+4(iy_i ) _ 1(Q + W-1,1,0 + wi ,-i,o + ™i,i,o + W-1,0,-1 + w -i,o,i + wi,o,-i + ^1,0,1 

+W ,-1 -1 + ^0 -1,1 + W ,l-l + U>0,l,l) 
+4 2 (w_i A0 + U>1,0,0 + ^0 -1,0 + Wo,l,0 + ^0,0,-1 + M> ,o,l) 

+4 3 w , ,o} • 

From the forms it is easy to conjecture the general form. 
If we set 

2 

W = F(xij X2, • • • , X n ) = ^ ^ Qiiia—tn^i ^"2 ' ' ' %n (^2) 



then 



ai+Aj /-a2+A 2 /-a n +A n 

/ ' ' ' / wdx\d%2 ■ ■ ■ dx n 

ai— Ai JCV2-A2 Join— A n 

= AlA2 3 :" A " £ )o^ 1J2 ..,„ (23) 

J1J2,- j»e{-i,o,i} 

where H{ji, j 2 , ■ • • , Jn}o is the number of in j 2 , • ■ • , jn} and 

w hh-u = F ( a i + J1A1, a 2 + j 2 A 2 , • ■ ■ , a n + j„A„) 

= /(«1 + JlAi, « 2 + J2A2, ■ • • ,a n + jrAn)- 

The formula is beautiful enough. 
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3 Concluding Remarks 

In the paper we gave a multidimensional analogue of the Simpson's formula of integral. As 
far as we know the formula has not been given, which is a bit mysterious. We believe that 
it will become useful in all fields related to integrals. 

Next let us consider more general problems. Let (M, g) be a curved space-time and / 
be a function on M. If we consider the integral 

f(x)dv g (x) 

M 

where dv g is some measure on M, it is almost impossible to calculate. For that we want to 
discretize the integral as follows. 

W /(,)**(,) ^ £ { some formula like {23} 

{D} {D} 

Therefore we must choose D to be "computable", which is just key point. 
Details of calculation will be reported in [2] . 
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